function __test_simpson()

	%Functie de test pentru metoda Simpson de calcul 
% a unei integrale pe intervalul [a b]
	%Testul nr.1
	I=simpson(-1, 1, 100, @fs1);
	printf("Valoarea integralei pentru functia exp(x) este :\n");
	disp(I);
        % solutia este: 2.3504
	
	%Testul nr.2
	I=simpson(-pi, pi, 100, @fs2);
	printf("Valoarea integralei pentru functia cos(x) este :\n");
	disp(I);
        % solutia este: 0
	
	%Testul nr.3
	I=simpson(-pi/4, pi/4, 100, @fs3);
	printf("Valoarea integralei pentru functia 1/(1+x^2) este :\n");
	disp(I);
        % solutia este: 1.3315

endfunction
	
